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ABSTRACT 

We study the prospects for observing H2 emission during the assembly of primordial 
molecular cloud kernels. The primordial molecular cloud cores, which resemble those at the 
present epoch, can emerge around 1 + z = 20 according to recent numerical simulations. 
The kernels form inside the cores, and the first stars will appear inside the kernels. A kernel 
typically contracts to form one of the first generation stars with an accretion rate that is as 
large as ~ O.O1M year -1 . This occurs due to the primordial abundances that result in a 
kernel temperature of order 1000K, and the collapsing kernel emits H2 line radiation at a 
rate ~ 10 35 erg sec -1 . Principally J = 5 — 3 (v=0) rotational emission of H 2 is expected. At 
redshift 1 + z — 20, the expected flux is ~ 0.01 /xJy for a single kernel. While an individual 
object is not observable by any facilities available in the near future, the expected assembly 
of primordial star clusters on sub-galactic scales can result in fluxes at the sub-mJy level. 
This is marginally observable with ASTRO-F. We also examine the rotational J = 2 — 
(v=0) and vibrational 5v — 1 emission lines. The former may possibly be detectable with 
ALMA. 

Key words: cosmology: observations — galaxies: formation — infrared: galaxies - 
ISM: molecules — submillimeter 
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1. INTRODUCTION 

Observation of the first generation of stars presents one of the most exciting challenges 
in astrophysics and cosmology. Hydrogen molecules (hLjs) play an important role as a cooling 
agent of gravitationally contracting primordial gas (Saslaw & Ziploy 1967) in the process of 
primordial star formation. It has previously been argued that H 2 is an effective coolant for 
the formation of a globular cluster (Peebles & Dicke 1969), if such objects precede galaxy 
formation. The contraction of primordial gas was also studied in the pioneering work of 
Matsuda, Sato, & Takeda (1969). The observational feedback from such Population III 
objects was examined in Carr, Bond, & Arnett (1984). In the context of the CDM (Cold 
Dark Matter) scenario for cosmic structure formation, Couchman & Rees (1986) suggested 
that the feedback from the first structure is not negligible for e.g. the reionization of the 
Universe and the Jeans mass at this epoch (e.g. Haiman, Thoul, & Loeb 1996; Gnedin & 
Ostriker 1997; Ferrara 1998). Recent progress on the role of primordial stars formation in 
structure formation is reviewed in Nishi et al. (1998). 

More recently, the first structures in the Universe have been studied by means of very 
high resolution numerical simulations (Abel, Bryan, & Norman 2000). The numerical resolu- 
tion is sufficient to study the formation of the first generation of molecular clouds. According 
to their results, a molecular cloud emerges with a mass of ~ 10 5 solar masses as a result 
of the merging of small clumps which trace the initial perturbations for cosmic structure 
formation. Due to the cooling of H2, a small and cold prestellar object appears inside the 
primordial molecular cloud. It resembles the core of molecular clouds at the present epoch. 
Their numerical results are consistent with other numerical simulations by Bromm, Coppi, 
& Larson (1999), whose results are consistent with those of Abel et al. (Bromm, Coppi, & 
Larson 2001). All the simulations predict that the primordial molecular clouds and their 
cores appear at the epoch of 1 + z ~ 20 (z is redshift). 

According to these results, the first generation of young stellar objects has a mass of ~ 
200 solar mass, density of ~ 10 5 cm -3 , and temperature of ~ 200 K. We stress that these are 
cloud cores, and not stars. The resolution of the existing numerical simulations is inadequate 
to extend the dynamic range to actual star formation. Resort must be made to analytical 
arguments. 

Inside the cores, a very dense structure appears. We shall call this a kernel for clarity 
of presentation. Its density increases to a value as high as ~ 10 8 cm~ 3 where the three-body 
reaction for H 2 formation occurs. Further evolution is by fragmentation (Palla, Salpeter, & 
Stahler 1983) and/or collapse (Omukai & Nishi 1998). In the current paper, we consider 
the case that the kernel collapses: the alternative case of fragmentation will be discussed 
in another paper. Recent theoretical work has suggested that fragmentation occurs if the 
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metallicity is above 10 3 5 x Z Q where Z & is the solar metallicity (Bromm et al. 2001). 

Omukai & Nishi (1998) show that the kernel collapses with a very large accretion rate. 
It is about 0.01 M year -1 . This accretion rate is so high due to the cooling of H 2 . The H 2 
rotational and vibrational level transitions have wavelengths in the the rest frame infrared 
band. The recent numerical simulations suggest that these are emitted at 1 + z ~ 20. The 
redshifted wavelength is detectable in the infrared - sub-mm bands. Of course, the predicted 
luminosity of individual kernels would not be expected to be extraordinary huge, although 
if the kernels are clustered, the assembly of kernels may be detected by future observational 
facilities. Hence, in this paper, we will consider observationally feasible conditions for the 
assembly of primordial kernels. Future telescope projects that are relevant are ALMA for 
the sub-mm range and ASTRO-F for the infrared range. 

In §2, we estimate the expected luminosity from a single contracting kernel. In §3, we 
confirmed that the luminosity from the kernel is primarily due to the emission of H2, and 
we examine which lines are dominant. In §4, the observational feasibility is discussed. 



2. LUMINOSITY OF PRIMORDIAL MOLECULAR CLOUD KERNEL 

According to the recent numerical simulations, molecular cloud cores appear prior to 
the formation of population III objects. Cores contain a dense and cool kernel, which has 
nearly the same physical properties as the initial conditions chosen by Omukai & Nishi 
(1998). Omukai & Nishi find that the kernel collapses in a freely falling manner and that 
the first prestellar object grows with a very high accretion rate. This is possible because of 
the cooling by line emission of H 2 . As a result, they find that the accretion rate is ~ 10~ 2 
M Q year -1 . 

We are interested here in the possibility observing primordial prestellar objects. We 
express the accretion rate in terms of luminosity. It is convenient for us to estimate the rate 
of energy release owing to the high primordial accretion rate. To determine this, we must 
examine the gravitational energy where the accretion occurs. To examine the gravitational 
potential energy, we need to determine the mass distribution around the centre of the ker- 
nel, where the first star emerges. Also, we assume a spherical configuration for the mass 
distribution. According to Omukai & Nishi, a high accretion rate is realized if a similarity 
collapse occurs with the adiabatic heat ratio of 1.1 (e.g. Suto & Silk 1988). The density 
distribution is described as 

gMr) -2 
<91nr 2-7' 1 ' 

Here, r is the radial distance from the centre, p(r) is the mass density of atomic H, H 2 and 
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He, and 7 is the specific heat ratio. We set the mass-density distribution of a protostellar- 
kernel with 7 = 1.1 as Pk( r ) = Pko( r / r o)~ 2 ' 2 where p k0 is 2.0xl0" 16 g cm" 3 and r is 0.01 pc. 
These values are appropriate for fitting the protostellar kernel of Omukai & Nishi. 

Once the density distribution of the kernel is given, we can estimate the gravitational 
potential at a radius r: $(r) = GM(r)/r where G is the gravitational constant and M(r) 
is the mass inside the radius r. Since the density is proportional to r -2 ' 2 , we find $(r) = 
GM(r)/r = 47rGpko r o' 2 /0.8r a2 . We are interested in the mass within r , then, $(r ) = 
2.0 x 10 11 erg mass -1 is obtained. Hence, as found from the dimension of "erg mass" 1 " for 
the potential, we obtain the energy release rate corresponding to the accretion rate as being 

L m = M x * (ro ) = 1.3 x 10-erg x (2) 
As found in the next section, this release of energy is explained by the cooling of H 2 radiation. 



3. H 2 EMISSION LINES 

The large accretion rate of the first stellar object arises because of the cooling of H 2 . 
Here, we also re-formulate the energy loss owing to line emission of H 2 . To estimate the 
cooling by H 2 , we should know the number fraction of H 2 , / 2 (r), and the temperature of gas, 
T(r), in the kernel. Both should have some spatial dependence. We use the fitting formula 
for Omukai (2000) which is adapted from Omukai & Nishi (1998). For f 2 (r) (solid line in 
Fig.l); 

Mr) = 0.0001 + 0.495 x V7V7 (3) 

e-X-P ^i n o cm -3 J I" exp ^ 10 n.o cm -3 J 

where n(r) = 10.0 8 cm" 3 (r/0.01 pc)" 2 2 (the maximum of / 2 (r) is set to be 0.5 by definition). 
For T(r) (dashed line in Fig.l); 



i 



T(r) = 1000 K ( —^1—] 15 . (4) 
v > V 10 cm " / 

According to Omukai (2000), molecular cooling can be dominant below a number density 
of ~ 10 15 '°cm" 3 . According to our estimates, the temperature is always smaller than 1000 
K. We usually set T 3 = T/1000 K. Line emission of H 2 occurs due to the changes among 
rotation and vibration states. As long as the large accretion occurs in the range 200 K - 
1000 K (Omukai 2000), we can safely adopt the formulation of Hollenbach & McKee (1979) 
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for rotational and vibrational emission of H 2 . Adopting their notation: 



9.5 x 10- 22 T 3 376 



—exp 



1 + 0.127: 



we estimate the cooling rate as 



0.13 



3.0 x 10" 24 exp 



0.51 



erg s , 



(5) 



A(rot) = n H2 7 r (l + CHr) 1 + n H2 7 r (l + Chv) 1 erg 



cm 3 s" 



(6) 



where Chi = n Hc d(rot)/n H , CH 2 r = ™H 2 cd(rot)/n H2 , n Hc d(rot) = Aj/rf, n H2C d(rot) = A//7j 2 , 
and A j is the Einstein A value for the J to J — 2 transition; 7j is the collisional de-excitation 
rate coefficient due to neutral hydrogen; and 7j 2 is that due to molecular hydrogen. The 
first term of L r denotes the cooling coefficient due to the higher rotation level (J > 2) and 
the second one due to J = 2 — >• transition. The vibrational levels of both terms are set to 
be zero. For our parameters of interest (250 < T < 1000; and n > 10 8 cm" 3 ), both (s are 
much smaller than unity for dominant levels of J because of the low temperature. This is 
because the maximum Aj values are at most 3.0 x 10~ 7 sec -1 . As a result, these correction 
factors are found to be 1 + ( ~ 1.0. 

For the vibrational transitions; 



6.7 x 10~ 19 exp 



then, we get cooling rate as being 



5.86 



+ 1.6 x 10~ 18 exp 



11.7 



erg s" 



(7) 



A(vib) = n R2 L v (l + Ch v ) 1 + ™h 2 7 v (1 + Ch 2 v) 1 erg cm 3 s 
^Hcd(vib)/n H , and Ch 2 v = n H2C d(vib)/n H2 . Here, n Hcd (vib) 
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where Ch v = n Hc d(vib)/n H , and Ch 2 v = n H2C d(vib)/n H2 . Here, n Hc d(vib) = Aj/^j, and 
nn 2 cd(vib) = where A,- L j is the Einstein A value for the i to j transition; 7^ is 

collisional de-excitation rate coefficient due to neutral hydrogen; and 7^ 2 is that due to 
molecular hydrogen. In our formula, only the levels of v — 0, 1 and 2 are considered. This 
is sufficient since the temperature is lower than 1000 K. The first term of L v is a cooling 
coefficient of Sv — 1 and the second term is that of Sv — 2. The second term has no effect 
on our conclusion because of the low temperature. Combining Eq.(6) and Eq.(8), we obtain 
the total cooling rate as A thin = A(rot) + A(vib) erg cm -3 s~ x in the optically thin regime. 

We need a cooling function which can be used in the optically thick regime. Then, 
adopting the following extension of 7 thin , we estimate the cooling rate as: 



■ thick 



^thin 1 -exp(-r) 



exp(-r cont ) 



(9) 
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where L thin is representative of each term of L r and L v , and r is determined for each L th " 
The effect of the continuum absorption below 2000 K is denoted by exp(— r cont ), in which 

-0.9 



Tcont = P(r)\j 



4.1 ( -L - - ) ^3" 4 ' 5 + 0.012p a51 (r)T q 2 - 5 
1 p(r) p 



(10) 



according to the estimate of Lenzuni, Chernoff, & Salpeter (1991) who obtain a fitting 
formula for the Rossland mean opacity in a zero-metallicity gas. Here, Aj is the Jeans length 
and p is 0.8 g cm -3 . Their fitting formula is reasonable if we consider the temperature 
range T > 1000 K. Our lowest temperature of the collapsing kernel is about 300 K. Then, 
their formula may not be appropriate, while it gives an upper limit if we adopt r cont of 1000 
K instead of really having r cont below 1000 K. For all of our estimates, r cont is sufficiently 
smaller than unity. Then, we can neglect the effect of continuum absorption. 

We shall estimate r for each line. Defining the suffix of % as the upper energy level and 
j as the lower energy level, we define as r for each line emission. 

Ajj(? n(x,i)gj Aj 
Tl3 8*1$. 9j 28v D 

Here, Aij is a transition rate from the % state to the j state; Uij is the central frequency 
of emission in the rest frame; c is the speed of light; ^ and gj are statistical weights for i 
and j respectively; Aj is the Jeans length ([2/c b T(r)//i(r)m H ] ' 5 x [3n /32Gp(r)]~ ' 5 ); k h is 
the Boltzmann constant; p(r) is mean molecular weight; and 8vd = [2/cb7 1 (r)//z(r)mH] ' 5 
is the velocity dispersion. The procedure of Jeans length shielding is useful for a simple 
analytical analysis (Low & Lynden-Bell 1976; Silk 1977). The population density of n(x,i) 
for rotational level is estimated with the statistical weight of gj = 2 J + 1 (i.e.i — > J). The 
suffix of x denotes each species. We consider n(H 2 , J) up to the upper level of J=5, which is 
selected since J = 5 is a dominant rotational level for our fitting formulae of Eq.(l), Eq.(3) 
and Eq.(4) as checked below. Then, we assume n(H 2 , J > 2)/n(H 2 , J = 2) = 27/5 

In our calculation to obtain Eq.(ll), for a rotational transition with v — 0, A 2t o = 2.94 x 
10~ n sec" 1 ; and Ajj_ 2 = 5A 2 , /162 x J(J - 1)(2J- l) 4 /(2j + 1) sec" 1 are considered. For a 
vibrational transition, A w = 8.3 x 10~ 7 sec -1 ; A 2 i — 1.1 x 10~ 6 sec -1 ; and A 2 o = 4.1 x 10 -7 
sec -1 are considered. 

To find the dominant rotational emission, we calculate 

^ = / lmnc i47rr 2 n(r)J(T(r))/ 2 (r)dr 
J kernr i47rr 2 n(r)f 2 (r)dr 



where 



7 T(r) x "■ r> 



J(T(r)) = ( - X ] . (13) 
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Here, J(T(r)) means the dominantly contributing J- level to the statistical weight at tem- 
perature of T(r) (Silk 1983). Using our fitting formula of /2(f) and T(r), we find J max is 
about 5. Then, we regard that the dominant rotational line cooling due to the first term of 
Eq.(5) is the J = 5 — 3 transition. To estimate the luminosity via L thlck of Eq.(9), hence, we 
use A$ : 3 since the J = 5 (v = 0) level is dominant. For J = 2 — (v = 0) transition (the 
second term of Eq.(5)), we use A 2t o to find L thlck . The second term of Eq.(7) is always much 
smaller than the other terms, and we ignore it here. 

The estimated total luminosity ( J kcrnel 4:iir 2 Adr) with modified L thlck s is 6.7 x 10 34 erg 
s _1 . This is comparable to the luminosity estimated in §2. Thus, we confirm that the large 
accretion rate is possible owing to H 2 cooling. The emission fraction is 0.82 for the rotational 
transitions with J > 2; 0.12 for rotational transitions of J = 2 — (v = 0); and 0.06 for 
vibrational transition of Sv — 1. Note that the fraction of J > 2 is not much larger than 
the fraction of J = 2 — (v = 0) transition. This is caused by the optical depth effect. If 
we adopted unity for the escape probability (i.e. any r much smaller than unity), we would 
find the result that there is dominance of J > 2 rotational cooling. For the same reason, the 
fraction of vibrational cooling can be smaller than that of J = 2 — (v — 0). 

The line broadening is estimated to be ~ 5vr> in the dimension of velocity. Adopting 
this, the luminosity per Hz (T=1000 K is assumed) is the following; The rotation emission 
of J = 5 — 3 (3.1 xlO 13 Hz; 9.7xl0~ 3 mm) is 1.9 xlO 26 erg Hz -1 , the rotation emission of 
J = 2 — (l.lxlO 13 Hz; 2.8xl0~ 2 mm) is 7.9 xlO 25 erg Hz -1 , and the vibration emission of 
v = 1 - (1.2xl0 14 Hz; 2.4xl0" 3 mm) is 3.4 xlO 24 erg Hz" 1 . 



In the current paper, we need the observational frequency and wavelength which are 
redshifted. For the three lines; we get 1.6 xlO 12 Hz; 1.9 xlO -1 mm for J = 5 — 3, 0.53 xlO 12 
Hz; 5.6 xlO -1 mm for J = 2 — 0, and 6.1xl0 12 Hz; 4.9xl0~ 2 mm for 5v = 1. The adopted 
1 + z is 20 according to the results of recent numerical simulations. 

To find the observed flux, we determine the distance to the object. Then, we calculate 
it numerically from the following standard formula: 



where H is the Hubble parameter and used 75 km sec^ 1 Mpc^ 1 , JIa is the cosmological 
constant parameter, and Qm is the density parameter. We consider the case £7a + = 1- 
Adopting D 19 , we obtain the observed fluxes of each of the lines. The results are summarized 



4. 



OBSERVATIONAL FEASIBILITY 




(14) 
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in table 1. Although the line-broadening is estimated from u 5vb/c, we also correct it with 
redshift effect (uq is the central frequency). For each of the parameter sets of (J1a,^m)) we 
obtain D 19 = 0.62 x 10 10 pc (Q A = 0, ft M = 1), D 19 = 1.00 x 10 10 pc (Q A = 0.7, Q M = 0.3), 
and D w = 4.58 x 10 10 pc (Q A = 0.9, Qm — 0.1). According to Table 1, the rotational line 
of J = 5 - 4 (v=0) is 0.07 fiJy; 0.05 fiJy for J = 2 - (v=0); and 0.001 fiJy for Sv = 1 if 
Qm = 1- Thus, we find that a single kernel is not observable by any available observational 
facilities in the near future. The expected flux is too small. 

However we note that if the kernels are assembled on a galactic scale, the assembly can 
be detected by ongoing projects. We shall estimate the number of kernels in a primordial 
galaxy. Firstly, we must consider the lifetime of a H 2 emission kernel. It may be 10 4 years 
for a 100 M star at the specified accretion rate of 0.01M Q year" 1 . A 10 6 M Q cloud is 
expected to form 1000 such stars at a plausible efficiency. Taking its lifetime to be 10 5 years 
as the free-fall time of a 1O 6 M cloud, we find that its luminosity is 100 L acc (see Eq.(2) for 
definition of L acc ). 

Next, we consider an entire primordial galaxy with lO n M . It may form 10 9 supernovae 
over its entire lifetime. We assume that it makes 10 7 primordial massive stars, since this 
number of massive stars gives enrichment to roughly 1 percent of solar metallicity. If the 
burst of formation of primordial molecular cloud kernels occurs in the central region of the 
primordial galaxy (we postulate 1 kpc as the size of the kernel-forming region), then there are 
10 4 such giant molecular clouds with 1O 6 M in the kernel-forming region. During this phase, 
the cumulative luminosity would be 10 6 L acc . It seems that the predicted line luminosity 
would be easy to observe. However, the dynamical time-scale in the kernel forming region 
may be 10 7 years (e.g. the duration of the starburst). Then, we obtain 10 4 L acc for the 
luminosity of the primordial galaxy undergoing its first star formation burst, conservatively 
assuming 10 5 years as the life-time of giant molecular clouds. Finally, for H 2 emitting proto- 
galaxies, we obtain 0.7 mJy for J = 5 — 3, 0.5 mJy for J = 2 — and 0.01 mJy for 5v — 1 
(Mm = 1.0). 

Thus, the rotational emissions of J = 5 — 3 and J = 2 — are at sub-mJy flux levels. 
Those frequencies enter into the sub-mm and far-infrared bands if the kernels form at 1 + z — 
20. These are observational ranges of ALMA and ASTRO-F. According to the current status 
of the instrumentation for ALMA, 80-890 GHz is the allowed detectable range. It is feasible 
for our prediction of the redshifted H 2 emission of the J = 2 — line. Unfortunately, however, 
the transmission is bad for the predicted feature at 530 GHz because of the atmosphere of 
the Earth. Then, it may be necessary to detect H 2 emission from a protogalaxy forming 
later than 1 + z = 20. We note that £ 1 + z — 20' is obtained by taking Vl M = 1- If < 1, 
which seems to be reasonable if Qa = 0.7, the formation of kernels is delayed. Then, we 
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can expect to detect emission at a higher frequency than 530 GHz and detection should be 
feasible of the H 2 emmision from the assembly of the primordial molecular kernels. 

For the case of the J = 5 — 3 line, this is observable at 194 /iin wave-length. This is 
in the allowed range of ASTRO-F (2-200 //m) . At the wavelength of 200 //m, unfortunately, 
it is difficult to detect sub mJy-level fluxes because of the properties of the detectors. If a 
higher rotational transition than J = 5 — 3 dominates, the observations are more suitable 
for the allowed range of ASTRO-F. The higher transitions emit more energetic photons, and 
enter into the possible range of ASTRO-F. This may be possible if the external radiation 
field (e.g. Susa & Umemura 2000) has positive feedback on the primordial kernels. Thus, 
although we need a very specific condition in order to detect H 2 emitting protogalaxies, it 
may be possible to achieve this goal with future observational facilities. 

Here, we advocate a deep blank field survey. The unit area for the survey can be 
estimated by the number density of QSOs, since QSOs are believed to form at high density 
peaks during cosmic structure formation. The burst of primordial kernel formation may be 
expected in the environment where QSOs can form. The most recent and reliable number 
density of QSOs is determined by Miyaji, Hasinger & Schmidt (2000). According to their 
luminosity function of soft X-ray AGNs, a typical number density (e.g. their figure 11) is 
about 10 -6 h\ Mpc -3 beyond z ~ 3. They have assumed that Hubble parameter is 50 
km/s/Mpc. Then, a unit covering area for a single QSO at a fixed redshift (i.e. a fixed 
distance) is typically 10 4 h^Q Mpc 2 . The proper transverse size of an area, -D S i ze , as seen by 
us, is its comoving size times the scale correction factor at the time of emission; 56D\g/20 
at 1 + z — 20, where 89 has the units of radians. Hence the unit area needed for our blank 
field survey is estimated to be 7.6 x 7.6 degree 2 if D sizc = 66.6 Mpc, D 19 = 10 10 pc, and 
H = 75. The field of view of ALMA is about 10 x 10 arcsec 2 . Thus, we need a huge 
amount of observing time to survey the expected area. The field of ASTRO-F for the long 
wavelength detector is 12.5 x 2.5 arcmin 2 . Although ASTRO-F plans to survey nearly all 
the sky, the detection limit for the survey mode is at 10 mJy level and we cannot expect to 
detect H 2 emission. We need good pointing and long exposure-time for the entire area of 
7.6 x 7.6 degree 2 . Then, the pointing mode of ASTRO-F may be suitable. It can detect line 
emission at the mJy level in a single pointing exposure. To detect the predicted sub-mJy 
level flux, we need to expose the same region several times. As a result, we can expect that 
an H 2 emitting protogalaxy could be found near the oldest QSO via deep and off-centred 
(i.e. avoiding saturation) spectroscopy using the pointing mode of ASTRO-F. 

It might be difficult to justify the necessary observing time for a blank field survey. In 
such a situation, we need to specify the observational target. Which object is suitable for 
our first observations? One possibility would be unidentified SCUBA sources (e.g. Smail, 
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Ivison, & Blain 1997). These are believed to have been detected at high redshift owing to 
the so-called "negative k-correction" . As long as they are primordial galaxies, the possible 
range of z for SCUBA objects is from 3 to 10. If z = 10-SCUBA objects are observed by 
near-future facilities, it may be possible to detect emission from H 2 at the sub-mJy level. 
But, in this case, the background radiation from dust prevents us from confirming H 2 line 
emission. To resolve this, a high dispersion spectrograph is needed. Furthermore, if the 
SCUBA objects are metal-rich primordial galaxies, the primordial molecular cloud evolution 
favours the fragmentation scenario (Bromm et al. 2001), and low mass star formation. 
The criterion of Bromm et al. for being at metallicity below 10~ 35 x Z Q to allow massive 
primordial star formation is only marginally satisfied at this redshift, as would be required 
to allow dominant H 2 cooling as suggested by Omukai (2000). 



5. SUMMARY 

One of the main goals of cosmology is to find the first generation of stars. When they 
form, strong H 2 emission is expected. We have examined the observational feasibility of 
detection. According to our analysis, J = 5 — 3 and J = 2 — rotational emission are 
expected (at v — 0). Protogalaxies can have sub-mJy fluxes when the first stars form at 
1 + z = 20 as an assembly. It seems that these systems are marginally detectable by ALMA 
and ASTRO-F. However, the transmissivity of air for the emission of J = 2 — is low 
for the ALMA project, and the wavelength edge of the filter for the far-infrared band is 
not sufficiently sensitive for the J = 5 — 3 line. Thus if future telescopes can detect H 2 
emission from primordial molecular cloud kernels, the kernels and first stars have formed 
after 1 + z = 20 and/or the temperature of the kernel should be larger than predicted in the 
simplest models. 
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Table 1: Expected Emission Lines 





J = 5-3 


J = 2-0 


5v = 1 






D 19 (10 10 pc) 


v (10 13 Hz; z = 0) 


3.09 


1.06 


12.21 








v (10 12 Hz; z = 19) 


1.55 


0.53 


6.11 








A (1(T 3 mm; z = 0) 


9.71 


28.29 


2.46 








A (1(T 3 mm; z = 19) 


194.2 


565.8 


49.2 








L u (10 26 erg sec" 1 Hz" 1 ) 


1.87 


0.79 


0.034 








/„ (10- 8 Jy; z = 19) ' 


7.34 


5.32 


0.148 


1.0 


0.0 


0.62 


/„ (10- 8 Jy; z = 19) 


2.85 


2.07 


0.058 


0.3 


0.7 


1.00 


/„ (10- 8 Jy; z = 19) 


1.29 


0.94 


0.026 


0.1 


0.9 


1.48 
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FIGURE CAPTIONS 

Fig.l - - Dcnsity-H 2 fraction relation (solid line) and Density-Temperature relation 
(dashed line) are depicted. Using Eq.(l), we translate them to Radius-Density and - 
Temperature relations, respectively. 



